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Abstract 

We consider cooling of neutron stars (NSs) with superfiuid cores composed of neutrons, 
protons, and electrons (assuming singlet-state pairing of protons, and triplet-state pairing of 
neutrons). We mainly focus on (nonstandard) triplet-state pairing of neutrons with the \mj\ = 2 
projection of the total angular momentum of Cooper pairs onto quantization axis. The specific 
feature of this pairing is that it leads to a power-law (nonexponential) reduction of the emissivity 
of the main neutrino processes by neutron superfluidity. For a wide range of neutron critical 
temperatures T cn , the cooling of NSs with the \mj\ = 2 superfluidity is either the same as the 
cooling with the m,j — superfluidity, considered in the majority of papers, or much faster. The 
cooling of NSs with density dependent critical temperatures T cn (p) and T cp {p) can be imitated 
by the cooling of the NSs with some effective critical temperatures T cn and T cp constant over NS 
cores. The hypothesis of strong neutron superfluidity with \mj\ = 2 is inconsistent with current 
observations of thermal emission from NSs, but the hypothesis of weak neutron superfluidity of 
any type does not contradict to observations. 

1 Introduction 



The cooling of neutron stars (NSs) strongly depends on the properties of the matter in their cores, 
primarily on the equation of state for the matter (Lattimer and Prakash 2001) and on the critical 
temperatures of NS nucleon superfluidity (Lombardo and Schulze 2001). At present, both are known 
incompletely, because there is no reliable microscopic theory. The available knowledge about the 
properties of the matter in NS cores can be improved by comparing observational data with NS 
cooling models. 

Here, we consider some cooling models. For definiteness, we assume the standard composition 
of NS cores: neutrons (n), protons (p), and electrons (e). 

It is well known that the neutrons and protons in NS cores can be in a superfiuid state (Yakovlev 
et al. 1999; Lombardo and Schulze 2001). Model calculations of superfiuid gaps show that proton 
pairing takes place in the singlet 1 So state of the proton pair. Neutron pairing can take place in 
the singlet and triplet ( 3 P2) states; singlet pairing generally arises in matter of moderate density 
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( p < po> where po = 2.8 x 10 14 g cm~ 3 is the nuclear density) and triplet pairing arises in denser 
matter. The triplet pairing, in turn, can be of several types that differ by the component of the 
total moment of the nucleon pair along the quantization axis (\mj\ = 0, 1, 2). Following the review 
articles by Yakovlev et al. (1999, 2001), we call the singlet pairing case A, the triplet pairing with 
\mj\ = case B, and the triplet pairing with \mj\ = 2 case C. Case C stands out among the 
remaining cases in that the superfluid gap in the dispersion relation for neutrons becomes zero at 
some points of the Fermi surface (at its poles). This leads to a fundamentally different (power- 
law rather than exponential) dependence of the NS neutrino energy losses on the gap amplitude 
(Yakovlev et al. 1999). The type of superfluidity with a minimum free energy occurs in nature. 
The A-type proton pairing in the NS core, the A-type neutron pairing in the crust and in the outer 
core, and the B-type pairing in the inner core are commonly considered in cooling calculations. 
The C-type neutron pairing in the NS core seems less realistic, but it is not ruled out by current 
microscopic theories. 

Up until now, virtually nobody has modeled the NS cooling with the C-type neutron super- 
fluidity. We know only the recent paper by Schaab et al. (1998), who attempted to take into 
account the effect of such pairing on the main processes of neutrino energy release in NS cores 
and on the NS cooling. Unfortunately, the authors failed to completely consider all of the factors 
that determine the thermal evolution of NSs. Thus, for example, when considering one of the most 
important neutrino processes, the neutrino emission during Cooper neutron pairing, the authors 
used an exponential dependence instead of the correct power-law dependence of the rate of energy 
release on T/T cn (T cn is the critical neutron temperature). 

Here, we analyze the NS cooling for the C-type neutron pairing more rigorously. We compare 
the coolings for the B- and C-type superfluidities. Our results are compared with observations. 

2 NS cooling for the C-type neutron pairing 

For the cooling to be adequately modeled, we must know the rate of neutrino energy release in 
superfluid NS interiors and the heat capacity of the NS matter. Recall that the main neutrino 
processes in NS cores are the direct Urea process, the modified Urea process, the nucleon-nucleon 
scattering reactions, and the Cooper pairing of nucleons. The first process, the most powerful 
mechanism of neutrino energy release, is a threshold process. The direct Urea process is open if 
PFn < 2ppp, where pp n and pf p are the neutron and proton Fermi momenta, respectively. 

The neutrino processes and heat capacity for cases A, B, and C have been studied extensively. 
The results can be found in the review article by Yakovlev et al. (2001) and in Gusakov (2002). 
Yakovlev et al. (2001) considered the neutrino processes with the A- or £>-type neutron and A-type 
proton superfluidities; they took into account the effect of the C-type superfluidity on the direct 
Urea process, the neutrino energy release during Cooper neutron pairing, and the neutron heat 
capacity. Gusakov (2002) investigated the modified Urea process during the C-type pairing and 
the nn- and np-scattering reactions. He also rigorously took into account the effect of the combined 
A- or B-type nucleon superfluidity on the modified Urea process. 

Here, we use the relativistic nonisothermal NS cooling code (see Yakovlev et al. 2001) adapted 
for the model with the C-type neutron superfluidity. This code allows us to construct the cooling 
curves, i.e., to determine the dependence of the surface temperature of a star, T£°, (recorded 
by a remote observer with allowance made for the gravitational redshift) on its age, t. Following 
Kaminker et al. (2002), we use the equation of state by Negele and Vautherin (1973) in the NS crust 
and the equation of state by Prakash et al. (1988) in the NS core with the compression modulus 
K = 240 MeV for symmetric nuclear matter at p = po (model I from their paper for the symmetry 
energy). For this equation of state, the direct Urea process opens at p > po = 7.851 x 10 14 
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g cm~ 3 , i.e., at M > 1.358Mq. In this case, the maximum NS mass is M = 1.977M Q . In the 
cooling code used, the effective proton and neutron masses in the inner NS layers renormalized with 
allowance made for multiparticle effects were assumed to be 0.7 of the mass of the free particles. 
The relationship between the NS surface and internal temperatures was taken from Potekhin et 
al. (1997). 

To better understand at which critical nucleon temperatures the C-type neutron superfluidity 
in NS cores differs from case B, we mentally replace the NS under consideration by a sphere with 
a radially constant density characteristic of this star. In that case, after the completion of thermal 
relaxation and before the beginning of the photon NS cooling stage, the time derivative of the 
internal NS temperature can be estimated as dT/dt ~ —Q/C, where Q is the energy released in 
the form of neutrinos per unit volume per unit time and C is the heat capacity per unit volume. 
This formula holds for both B- and C-type neutron superfluidities. Since the C-type superfluidity 
suppresses the NS energy release more weakly and since the energy released via Cooper pairing is 
higher than that in case B, we may say that Qc > Qb always. Similarly, Cq > Cb, therefore, it 
is not clear in advance which stars cool down faster: with the C- or U-type neutron superfluidity 
in the NS core. If we introduce the critical nucleon temperatures T cn and T cp for the characteristic 
density under consideration, then, depending on T cn and T cp , the function \og[QcC 'b / '(Q bCc)] will 
be convenient for understanding when the cooling for cases C and B proceeds almost identically 
and when the differences are at a maximum. 

Let us examine Fig. 1. It consists of six pictures, each corresponding to a certain temperature 
(T = 10 8 , 3 x 10 s , or 10 9 K) and density (p = 2po, or 5po). The direct Urea process is forbidden 
at p = 2po and permitted at p = 5po- Our calculations show that when the temperatures and 
densities are varied, the pictures do not change fundamentally until the threshold density po at 
which the direct Urea process switches on (or off) is crossed. The logarithms of the critical proton 
and neutron temperatures at a fixed density are plotted along the vertical and horizontal axes of 
each picture, respectively. The critical nucleon temperatures lower than 10 7 K (weak superfluidity) 
have virtually no effect on the NS cooling and, hence, are not considered. In the figure, the function 
^°s[QcCb /(QbCc)] is indicated by different shades of gray. The shade of gray for the region of 
critical temperatures to the left from the vertical dashed line (nonsuperfluid neutrons) corresponds 
to zero of this function. The darker is the shade of gray, the larger is log[QcCs / (Q bCc)]- F° r 
example, for T = 10 8 K, the maximum of this function reaches (6 - 7) for critical proton and 
neutron temperatures T cPiCn ~ (10 9 - 10 10 ) K at p = 2p and (8 - 9) for T cp ~ (10 7 - 10 8 ) K, 
T cn ~ 3 x (10 9 — 10 10 ) K at p = 5po. Conversely, the lighter is the shade of gray relative to the 
'zero' shade, the smaller is the function. The region with the shade of gray lighter than the 'zero' 
shade is narrow and lies near T cn ~ 3T 1 . In this region, the minimum value of the function is of 
the order of -1 for both densities. A further increase in the internal NS temperature causes the 
differences between the B- and C-type neutron superfluidities to be smoothed out; the smoothing 
occurs faster for the densities at which the direct Urea process is open. Thus, for T = 10 9 K, 
the maximum of loglQcCs / (Q bCc)] is (3 - 3.5) for p = 2po and (1 - 1.2) for p = 5po, while its 
minimum is nearly zero for both densities. Summing up our results, we may say that the inequality 
(dT/dt) c > (dT/dt) B holds for any T, T cn , T cp , and p. Therefore, one can hardly expect the 
NS cooling curve for case C to pass well above the cooling curve for case B (at fixed T cn and 
T cp ). Clearly, the reverse can be easily realized by an appropriate (see above) choice of critical 
temperatures. 

Let us now turn directly to an examination of the NS cooling curves. For simplicity, the neutrons 
in the stellar crust are assumed to be nonsuperfluid. This assumption has no effect on the difference 
between the cooling curves for the B- and C-type neutron superfluidities in the NS core. Although 
the critical nucleon temperatures T cn (p) and T cp (p) are actually functions of the density and vary 




Figure 1: The function log [QcC b / (Q bCc)} over a wide range of T cn and T cp is indicated by different 
shades of gray. The larger is the value of the function, the darker is the shade (see the text). 
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Figure 2: The cooling curves constructed for the B- and C-type neutron superfluidities in NS cores 
for various critical nucleon temperatures. The curves with identical T cn and T cp but for different 
superfluidities, B or C, are represented by the same type of lines (solid, dotted, and dashed lines). 
The long dashes correspond to the cooling without superfluidity. The NS mass is M = 1.35M Q . 



Table 1: Chosen nucleon superfluidity models. 



Number 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


11 


12 


W10 9 , K 


2.5 


5.0 


10.0 


3.2 


0.2 


0.5 


7.0 


0.85 


9.5 


0.2 


0.1 


5.0 


T cp /10 9 , K 


2.5 


0.25 


10.0 


0.6 


5.0 


1.0 


0.01 


3.0 


0.85 


0.7 


5.0 


0.01 



along the stellar radius, it was noticed that following certain semiempirical recipes, constant (over 
the core) critical temperatures can be matched to them. In this case, the cooling curves remain 
the same as if we took into account the exact dependence T cn (p) and T cp (p) (e.g., for protons, 
a constant effective temperature T cp close to T cp (p c ), where p c is the central density of the star, 
always corresponds to the temperature profile T cp (p)). Therefore, for simplicity, we assume the 
critical nucleon temperatures to be density-independent. 

Figures 2-5 show the cooling curves constructed for two NS masses: M = 1.35M (Figs. 2 and 
3) and M = 1.44Mq (Figs. 4 and 5). The direct Urea process is forbidden for M = 1.35M Q and 
permitted for M = 1.44Mq. In Table 1 we give the various critical nucleon temperatures T cn and 
T cp for which the cooling curves were constructed. Each pair of critical temperatures in the table 
has its own number. In Figs. 2-6, the cooling curves are marked by these numbers. The letter B 
or C near the number specifies the type of neutron superfluidity. For example, 3C means that the 
cooling curve was constructed for the C-type neutron superfluidity for the critical temperatures 
T cn and T cp under number 3. We see from analysis of Figs. 2-5 that all our main conclusions 
drawn from an examination of Fig. 1 are confirmed. For M = 1.35Mq, the largest differences 
are observed for the 3B and 3C curves in Fig. 2 (high critical temperatures). At critical proton 
temperatures T cp ~ 5 x 10 8 K, the neutrino emission during Cooper proton pairing is the most 
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powerful process in the NS core. Accordingly, the differences between the cooling curves for cases B 
and C decrease. This can be verified by comparing the AB and 4C curves in Fig. 3. As the critical 
proton temperatures decrease to T cp < 3 x 10 8 K and for T cn > 10 9 K (but T cn < 10 10 K), the 
differences again increase (the 2B and 2C curves in Fig. 2). When T cn > 10 10 K and increases, the 
differences gradually decrease and the cooling curves for the B- and C-type neutron paring coincide 
for T cn ~ 10 11 K. At moderately high critical temperatures T cn < 10 9 K, the curves for case B 
almost coincide with those for case C (see, e.g., Fig. 3: 5B and 5C or 6B and 6C). The kink in 
the 5B and 5C curves occurs at T = T cn and stems from the fact that the neutrino emission due to 
Cooper neutron pairing switches on. Note that the critical temperatures of the 6B and 6C curves 
exactly fall within the region where the temperature for the -B-type neutron superfluidity decreases 
faster than that for case C (log[C;cCB / '(QbCc)} < 0). This actually leads to the fact that the 
6B curve initially passes slightly below the 6C curve. Subsequently, as the internal temperature T 
decreases, the situation is reversed, because we reach the regions with faster cooling for the C-type 
superfluidity (see Fig. 1). 

For M = 1.44M©, the differences are largest for the 7B and 7C curves in Fig. 4, i.e., at 
high critical neutron temperatures and at low critical proton temperatures. As T cn decreases, the 
differences between the cooling curves for the B- and C-type superfluidities gradually disappear 
(see Fig. 4: 8B and 8C; Fig. 5: 10B and 10C). The remaining remarks almost completely repeat 
those for M = 1.35M . 

Our results can be of use in explaining observational data. The NS cooling with weak B- 
type neutron superfluidity was considered by Kaminker et al. (2002). The authors showed that the 
available observational data on the thermal radiation from eight isolated NSs could be explained by 
assuming strong proton (T cPjmax > 5 x 10 9 K) and weak neutron (T cnjmax < 10 s K) superfluidities. 
In this case, the critical temperatures T cp and T cn must depend on density p [have a maximum 
at p ~ (2 — 3)/9o> as predicted by microscopic theories of superfluidity; see, e.g., Lombardo and 
Schulze (2001)]. The mechanisms of strong proton and weak neutron superfluidities are primarily 
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needed to interpret the observations of two sources, RX J0822-43 (Zavlin et al. 1999) and PSR 
1055-52 (Ogelman 1995), the hottest ones for their age. The observational data on these sources 
are shown in Fig. 6 [we took the same T e °° and t as those used by Kaminker et al. (2002)]. Until 
recently, in addition to these two sources, yet another source, RX J185635-3754, must have been 
considered (Pons et al. 2002; Burwitz et al. 2001; Gansicke et al. 2001; Kaplan et al. 2001). 
However, in the just published paper by Walter and Lattimer (2002), the age of RX J185635-3754 
was revised (reduced by almost half) and is now estimated to be t ~ 5 x 10 5 yr, which is attributable 
to a decrease in the distance to this source. 

When comparing observations with cooling models, we should take into account the fact that 
the surface temperatures and ages of isolated cooling NSs are determined with a low accuracy [for 
the reasons discussed, e.g., by Pavlov and Zavlin (1998) and Yakovlev et al. (1999)]. The NS age 
t is probably known to within a factor ~ 2. 

Since the determination of T e °° and t is ambiguous, we can offer an interpretation of the observa- 
tions alternative to that offered by Kaminker et al. (2002). More specifically, if we take the lowest 
possible values of T e °° and t for PSR 1055-52, then the observations of the eight NSs considered 
by Kaminker et al. (2002) can be interpreted by assuming strong neutron superfluidity and weak 
proton superfluidity in the stellar cores. Although this interpretation is less plausible than that 
offered by Kaminker et al. (2002) (according to microscopic calculations of T cn {p) and T cp (p), the 
triplet neutron pairing is generally weaker than the single proton pairing) , it cannot be ruled out in 
advance and should be studied. It actually requires lower values of T e °° and t of the fairly old and 
hot source PSR 1055-52, because strong neutron superfluidity significantly reduces the NS heat 
capacity and speeds up the NS cooling at the late photon cooling stage [implying that the models 
by Kaminker et al. (2002) can describe older and hotter sources]. 

A detailed discussion of the hypothesis of strong neutron superfluidity and weak proton superflu- 
idity in NS cores is beyond the scope of this paper and will be given in our subsequent publication. 
Here, we only illustrate the possibility of this interpretation of observations for the two hottest 
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sources (RX J0822-43 and PSR 1055-52) by using simplified models of cooling NSs with critical 
temperatures T cn (p) and T cp (p) constant over the stellar core. The dotted lines in Fig. 6 represent 
the cooling curves for T cn = 5 x 10 9 K and T cp = 10 K. For comparison, the dashed line was drawn 
for strong proton superfluidity and weak neutron superfluidity (T cp = 5 x 10 9 K, T cn = 10 8 K). It 
virtually coincides with the cooling curve obtained by Kaminker et al. (2002) to interpret the above 
two sources. We see that the upper dotted line (S-type neutron superfluidity) almost coincide with 
the dashed line for t < 10 5 yr and passes slightly below the dashed line for t > 10 5 yr (because 
of the above suppression of the neutron heat capacity by strong neutron superfluidity) . Neverthe- 
less, given the inaccurate determination of T e °° and t, the upper dotted line could be assumed to 
satisfactorily explain these two sources. At the same time, the lower dotted line for strong C-type 
neutron superfluidity passes well below (because of the weaker suppression of neutrino emission by 
the C-type superfluidity) and cannot be considered to be acceptable. Therefore, the observations 
can (with the above reservations) be explained by strong B-type neutron superfluidity but cannot 
be explained by the C-type neutron superfluidity with the same T cn = 5 x 10 9 K. With a further 
increase in T cn , the upper dotted line would not change at all (strong B-type neutron superfluidity 
with T cn ~ 5 x 10 9 K almost completely switched off the neutrino reactions involving neutrons and 
made the neutron heat capacity equal to zero), while the lower dotted line would slowly approach 
the upper curve (the C-type neutron superfluidity increasingly suppresses the neutrino luminosity 
and heat capacity), reaching it at T cn ~ 10 11 K. Thus, the interpretation of the observations in the 
model of strong C-type neutron superfluidity is possible in principle but is implausible, because it 
requires unrealistically high values of T cn ~ 10 11 K. 

3 Conclusions 

The cooling of NSs with C-type neutron superfluidity in their cores has been consistently studied 
for the first time. The critical nucleon temperatures were assumed to be constant over the stellar 
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Figure 6: The cooling curves for the models with the B- and C-type neutron superfluidities in 
the core of a NS with M = 1.35M0 and observational data on the surface temperatures of two 
sources. The dashed curves are for T cn = 10 8 K and T cp = 5 x 10 9 K; the dotted curves are for 
T cn = 5 x 10 9 K and T cv = 10 7 K. 



core. A comparison with B-type neutron superfluidity was made. We found that the cooling curves 
for case C can pass well below the cooling curves for case B (faster cooling). Otherwise, the curves 
differ only slightly. The cooling with variable (over the core) critical temperatures T cn {p) and T cp {p) 
can generally be described by the cooling with some effective constant temperatures T cn and T cp . 

The model with C-type neutron superfluidity can account for the observational data on the 
thermal radiation from isolated NSs by assuming strong proton superfluidity (T cp > 5 x 10 9 K) 
and weak neutron superfluidity (T cn < 10 8 K). The interpretation of the observations in terms of 
the model with strong neutron superfluidity and with weak proton superfluidity (T cp < 10 8 K) is 
unlikely (in contrast to the same model for case B), because it requires too high critical neutron 
temperatures, T cn ~ 10 11 K. Our results can be useful in understanding which type of superfluidity 
(B or C) occurs in nature. 
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